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Fiber guided optical signal propagating in a Erbium doped nonlinear resonant medium is known to produce 
cleaner solitonic pulse, described by the self induced transparency (SIT) coupled to nonlinear Schrodinger 
equation. We discover two new possibilities hidden in its integrable structure, for amplification and control of 
the optical pulse. Using the variable soliton width permitted by the integrability of this model, the broadening 
pulse can be regulated by adjusting the initial population inversion of the dopant atoms. The effect can 
be enhanced by another innovative application of its constrained integrable hierarchy, proposing a system of 
multiple SIT media. These theoretical predictions are workable analytically in details, correcting a well known 
result. 

(M ■ I. INTRODUCTION 



Optical communication through fiber has achieved phenomenal development over the last two 



"^J- ', decades OptCom. Dissipation and dispersion in the media, which are the main hindrances in signal 



transmission, are usually attempted to be solved by the dispersion management techniques and de- 
vices OptConijdismang. On the other hand, in soliton based optical communication, mediated by the 



nonlinear Schrodinger (NLS) equation proposed much earlier nlssplit, the group velocity dispersion 
can be countered by the self phase modulation in the nonlinear fiber medium agarwal. However, the 
experiments revealed insufficiency of the model for its efficient practical application nlsdrab. Another 
proposal with improved solitonic transmission was due to the self-induced transparency (SIT), pro- 
duced by the coherent response of the medium to an ultra short optical pulse SIT, lamb. Finally, 
the benefits of both the NLS and the SIT systems, were combined in a coupled NLS-SIT model 
Maimitsov, by transmitting the optical soliton through an Erbium doped nonlinear resonant medium 
nakazawal makazawa2 . 

However the solitonic communication, in spite of its favorable features and theoretical advantages 
due to its underlying integrability, did not receive the needed response. Our aim here therefore, is 
to revisit the NLS-SIT model for exploring new possibilities hidden in its integrable structures and 
use them for the control and amplification of the solitonic pulse. Though the soliton usually moves 
with a constant velocity or speed, the integrable property of this coupled NLS-SIT system, as we find 
here, allows the soliton speed to be a tunable function. And since the soliton width in this case is 
related to its speed, which in turn is linked here to the initial population inversion, the pulse width can 
be regulated by manipulating the population inversion profile of the dopant atoms. This controlling 
effect can be further enhanced by exploring another specialty of this integrable system, namely its 
constrained integrable hierarchy, through a novel use of coupled multiple SIT media, in place of the 
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conventional single doped medium. The details can be worked out exactly due to the integrability of 
the model, detecting the limitation of a well known result on NLS-SIT soliton nakazawal,nakazawa2. 

Propagation of a stable optical pulse through a fiber medium, serving as a dispersive and nonlin- 
ear wave guide with Kerr nonlinearity nlssplit, can be described by the optical electric field E(z,t) 
satisfying the well known NLS equation 

iE z - E tt - 2\E\ 2 E = 0, (1) 

with space and time variables being interchanged as customary in nonlinear optics agarwal. On the 
other hand, ultra short optical pulse producing SIT in the medium can be described by the Maxwell- 
Bloch equation SIT, lamb 

%E Z = 2p, ip t = 2NE, iN t = ~(Ep* - E*p), (2) 

with the induced polarization p and the population inversion N of the medium, contributed by the 
Bloch equation. Fascinatingly, it is possible to combine these two effects by transmitting the stable 
nonlinear pulse produced in the fiber wave guide through a doped medium with coherent response, 
governed by a coupled NLS-SIT system given by a deformed NLS equation 

iE z -E tt -2\E\ 2 E = 2 Pl , (3) 

together with the SIT equations 

ipu = 2(N 1 E - won), iN lt = -(Epl - E* Pl ), (4) 

representing a nonholonomic constraint kundu09. In , p\ = vv* is the induced polarization and 
Aq = |z/| 2 — |z/| 2 , — 1 < N\ < 1, is the population inversion of the two level dopant atoms with nor- 
malized wave functions v, v for the ground and the excited states, respectively and wq is the natural 
frequency of these resonant ions. Assuming a homogeneous broadening of the frequency spread with 
a sharp resonance at Aw = w — wo, we have taken the symmetric distribution g(Aw) = 5 (Aw) and 
replaced the average value < p\ >= / dwg(Aw)pi(z,t,w) appearing in ([3|) by p\ = pi(z,t,wo) and 
normalized the coupling constants to ensure the integrability of the model. 

II. INTEGRABILITY AND SOLITON SOLUTION 

Recall that, an integrable nonlinear equation may be associated with a linear system <&t = 
J7$, $2 = V<&, defined through a Lax pair U(X), V(X), which are matrices with their elements depend- 
ing on the basic fields and a parameter A, known as the spectral parameter. Through compatibility of 
the linear Lax equations, inducing flatness condition: U z — Vt + [U, V] = 0, the Lax pair yield the given 
nonlinear equation and at the same time can be used to extract its exact solutions through the inverse 
scattering method (ISM) soliton. We remind again that, the space z and time t are interchanged here 
in the context of nonlinear optics. It is noteworthy that, while the set of coupled NLS-SIT equations 
([3]|1]) generalize both the NLS dH) and the SIT ([2]) equations, the associated Lax pair U n i s - S u, V n is-.sit 
contain these subsystems as the constituent parts: 

Uni S : S u(X) = U nls (\) = U S it(X), V n i s:sit (X) = V nls (X) + V si t(X) (5) 
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where U n i s , V n i s and U s u, V s u are the Lax pairs related to the NLS and the SIT equations, respectively. 
The NLS Lax pair is well known as sohton 

U nls (\) = i(a 3 X + U^), U<® =E<x + + E*a- (6) 
V nls (X) = V^ + V^X + V^X 2 , 

V® = (U^ - i(U<M)*y, y« = 2iU^\ V® = 2m 3 , (7) 

where a 3 , cj ± = \{o~ l ± ic 2 ) are the 2x2 Pauli matrices. The SIT Lax pair can be given by the same 
time-Lax operator as that of the NLS: U S it(X) = U n i s (X), while the space-Lax operator 

V s u(X) = i(X-w y 1 G 1 , G x = iVia 3 + Pl a + + p\a~ , (8) 

can be linked to the nonholonomic deformation kundu09. We check easily that, the flatness condition 
of (|6|7p yields the NLS equation ([I]), while (|6l8p the SIT equation ([2]) and similarly, the Lax pair ([5]) 
would yield the coupled NLS-SIT equations (|3|4j) 

The time-Lax operator U (A) plays the central role in the ISM for finding the exact solutions of the 
nonlinear equation soliton and therefore, since U(X) is the same for NLS (HJ, SIT ([2]) and NLS-SIT 
(|3l4p equations as seen from ([S]), the form of soliton solutions and the ISM procedure are remarkably 
similar for all the three equations. We therefore present soliton solutions for all of them in an unified 
way following the ISM sohton, which though an involved method, gives the 1-soliton solution in an 
amazingly simple form: 

E = c T^Tl2> 5 = exp[2i(u oo t + v oo z + 0)], (9) 

with c,(f) = constants. Note that, the crucial elements u 00 ,^ 00 = - J z dzv°° in Q, though linked to 
the Lax pair of the given system, need information only about their asymptotic properties: u°° = 
c 3 j7(Ai)|t=_ 00 , v°° = <7 3 V (Ai)|t = _oo, at discrete spectral parameter Ai. Therefore, fixing the initial 
condition of the basic fields involved in the NLS-SIT equations as 

E(z,t = -oo) -> 0, pi(z,t = -oo) 0,iVi(z,t = -oo) = N (z) (10) 

with an arbitrary function Nq(z), we can easily derive from the Lax pair (J5KHJ) : 

a3n nls:sit = U n i s:sit (Xi)\t=-oo = U n i s (Xi)\ t =-oo = U s it(Xi)\ t =-oo = 2i<7 3 Ai, 

^nis: Sl t = ^ :S it(A)| t= -oo = (^(A) + K/ s (A))| i= -oo = ^^ 3 (2A2 + (Ai- W ;o)" 1 A r o(^)), (11) 

which at the discrete spectral parameter with complex value: X± = k + ir/, take the explicit form 

<i s -,it = iAi, v% s .. sit = v^ s + v~ 4 , v% s = 2iXl v~ = i(Ai - wo^Noiz). (12) 

Inserting the needed complex valued expressions ()12p in (J9j> and grouping its real (Re) and imaginary 
(Im) parts we get the 1-soliton solution for the optical field in NLS-SIT equations (|3l4p in the familiar 
sech- form 

E = -lit] sech2Ce 2ie , ( = r](t - t - vz), 6 = ujz + kt + O , (13) 

where to an d 4>o are constant time and phase shift. Inverse speed v and phase rotation uj for the 
NLS-SIT soliton (fT3"P are given by the superposition 

V = V n l s + V sit , UJ = UJ n l s + UJ sit , (14) 
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of the corresponding parameters from the NLS and the SIT subsystems, derived from (|12p using 
Ai = k + ir\ as 

v nls = -~Im[v% s ] = -4k, u nls = Re[v% s ] = 2{k 2 - rf), (15) 
V 

v slt = --/m[v s °° t ] = -f{z), u sit = Re[v%] = -f(z), (16) 
VP P 

with k = k-w , p = k 2 + r] 2 , and f(z) = \f z N (z')dz'. 

It is intriguing to note that, since the z-evolution of the optical field E in the NLS-SIT model 
follows the superposition rule Q14|) contributed separately by the NLS and the SIT parts, the term iE z 
in equation (|5J), evolving according to solution (JT3J) , breaks up into two parts: one follows the NLS 
contribution with parameters (|15p and satisfies the pure NLS part of the equation in the left hand 
side, while the other part equates to the SIT deformation 2p± in the right hand side of ([3]) involving 
the related parameters (|16|) . Using this dynamics we derive the soliton solution for the dipole p\ from 
(|13|) . in the form 

p 1= 1 iV sech2C(^-r ? tanh2C)e 2 ^, (17) 
P 

with Q,9 as expressed in (fl"3|) . Inserting solutions (|13ll7j) for E and pi in (jlj) and integrating by t we 
derive further the solution for population inversion 

2 

Ni = N (l — — sech 2 20, (18) 

again in the solitonic form with arbitrary function Nq(z) = N\(t — > — oo), adjusted by the integration 
constant. We obtain thus the complete set of exact soliton solutions to the NLS-SIT equations ([3]- 
H]) as (fl"3|) for the optical field E, (fT7|) for the dipole pi and (fTKj) for the population inversion Aq. 
A beautiful interaction pattern can be noticed in these solutions, manifested in the superposition 
relations: v = v n i s +v s it, co = uj n i s +UJ s it, for the solitonic parametersappearing in (|13|14p . Intriguingly, 
in the absence of the SIT system with p\ = Aq = 0, when the coupled NLS-SIT equations reduce to 
the NLS equation ([T]) for the field E, one recovers from f)13f) the well known NLS soliton by simply 
putting v s u = oJsit = due to the vanishing of (jHJ). Therefore, the NLS soliton takes exactly the same 
form as (fT3j) . though the parameters are reduced to pure NLS case: v = v n i s , uj = ui n i s . Similarly, 
we can directly get the soliton solution for the pure SIT equations ([2]) in the same form (|13|17I18|) . 
but with soliton parameters reducing to v = v s u, uj = ui s it, due to switching off the NLS influence: 
Vnls = ^nis = 0- Thus our exact NLS-SIT soliton can reproduce the solutions for both the NLS and 
the SIT equations in a unified way, consistent with the ISM. However this rich interaction picture 
seems to have been missed in a well known earlier work nakazawal,nakazawa2, leading to wrong 
conclusions in the general case. In particular, the soliton solution for the NLS-SIT equation presented 
in nakazawal i nakazawa2 gives the expression for the pulse delay as 5 = — (1 + 7) ((4.9) in nakazawa2), 
which is equivalent to the inverse soliton speed for the SIT ((2.22) in nakazawa2), i.e. v = 5 = v s a, 
in our notation solNLSsit. Similarly, the phase rotation in nakazawal,nakazawa2 is given as a = 2n 2 , 
meaning u) = a = —ui n i s , (at k = 0, see (|15p ) in our notation solNLSsit. Both these results for the 
coupled NLS-SIT equations appear to be incomplete, when compared with our exact result (|14p . It 
is clear that, the solution of nakazawal,nakazawa2 can be justified only in a very limited sense, when 
v n is = uisit = and therefore unlike our soliton solution can not interpolate between the solutions of 
the NLS and the SIT equations. 
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This partial result unfortunately led to wrong conclusions, for the NLS-SIT system in general, stat- 
ing that (sect. IV nakazawa2), the normalized speed (i.e. of the NLS-SIT soliton is determined 
only by the SIT effect and similarly, the z dependence of the phase of the dipole (i.e. a) is determined 
solely by the nonlinear phase change due to the NLS soliton. Our exact solutions for the optical field E 
(fT3j) and the dipole p\ (JTTj) with correct expressions (fl^]) . conclude on the other hand that, only a part 
(i.e. v s it) in the normalized speed v^ 1 = {v n i s +v s i t )~ l of the NLS-SIT soliton is determined by the SIT 
effect, while there is an additional contribution coming from the NLS part v n i s . Similarly, the z depen- 
dence of the phase of the dipole and the input optical field gets contribution from both the NLS and 
the SIT parts as uj = uj n i s +oj s it, consistent with the interaction picture in the coupled NLS-SIT system. 

III. CONTROLLING OPTICAL SOLITON EXPLOITING INTEGRABLE STRUC- 
TURES 

Based on the integrable structures underlying the NLS-SIT system describing the propagation 
of optical soliton in fiber guided doped medium, we propose two possible ways for controlling the 
amplitude and width of the optical pulses. 

A. Soliton control by regulating initial population inversion 

It is commonly believed that, the exact soliton solution of a homogeneous equation always moves with 
a constant speed, width and frequency, as in the case of the NLS soliton (|15f) with constant values for 
v n i s ,u n i s . However, it is crucial to note that, for the NLS-SIT soliton the parameters (v,uj), as evident 
from (|14ll6p can become variable functions, depending on the initial population inversion Nq(z) (|10p . 
Due to this peculiarity of integrable structure of the NLS-SIT system, hidden in the expressions like 
(|5l8ll2ll6p . the soliton speed: v~ l and width: {rjv)~ l , as defined from the soliton argument Q (fT3|) . can 
be variable and linked to a controllable arbitrary function Nq(z). 

We show that, this important observation embedded in the integrability of the NLS-SIT system 
can open up a new avenue for controlling the optical soliton propagating through the doped medium, 
by regulating its initial population inversion profile Nq(z). This fact however remained unexplored 
in earlier investigations nakazawal,nakazawa2,kakei,porsezian, due to the restriction to a fixed initial 
profile Nq(z) = — 1. Note that, at this particular value giving f{z) = —1, our more general result (1161) 
reduces to the simplified expressions obtained earlier: 

lk 

Vsit = — i u sit = --, (19) 
P P 

The choice for the initial population inversion in the NLS-SIT model as an arbitrary function 
Nq(z) > —1, that we propose here, gives us the needed freedom for obtaining the excited and the 
ground state occupancies at the initial moment as |z>| 2 = ^(1 + Nq(z)) and \v\ 2 = |(1 — Nq(z)), 
respectively. Therefore, for Nq > —1, giving \D\ 2 > 0, we can prepare the dopant atoms initially in 
an excited state by optical prepumping, resulting to the creation of a laser-active amplifying medium 
with its intensity determined by Nq. Note that, only in such a case when more active dopant atoms 
are in the excited state, the optical soliton can gain net energy OptCom. 

In addition to the soliton pulse amplification, variable initial profile Nq(z) > —1, permitted by the 
integrability of the NLS-SIT system, can play a crucial role in controlling the shape and dynamics of the 
optical soliton. It is possible, as we see below, to address the important problem of pulse broadening 
by regulating the initial profile of the dopant atoms. For example, a solitonic pulse governed by the 
NLS equation under small perturbation by a term —i^E with T « 1, would suffer broadening 



5 



by a factor (4kr/(z))~ 1 , which can be worked out through the variational perturbation method as 
r)(z) = T] e _r ^agarwal, which is valid however upto the range Tz ~ 1 along z. Beyond this range with 
z » 1, as shown by some other method, the broadening of the pulse width follows a different rule, 
by increasing linearly with z at a rate slower than the linear medium agarwal. Though an attenuation 
with intensity loss would also occur simultaneously, the broadening leads to more serious problem of 
information loss and bandwidth limitation. Therefore we concentrate here only on the broadening 
problem of the perturbed NLS soliton, due to the increasing solitonic width e Tz along z, as shown 
in Fig 1. As stated above for z » 1 it would follow a different rule. We show that, by transmitting this 
solitonic pulse through a doped resonant medium, described by an interacting NLS-SIT model (|3]|1|), it 
is possible to control the pulse broadening, by suitably preparing the initial population inversion profile 
Nq(z). Fig 2a shows this controlling effect, where the broadening of the solitonic pulse suffered in Fig 
1, is countered by the narrowing of the pulse due to variable width V(z) = {y n i s + Nq(z) / (pT))^ 1 , by 
taking Nq(z) ~ ri(z)~ l . Note that the profile Nq(z) has to be adjusted differently at different ranges, 
as mentioned above, to control the broadening in the respective regions for a wide range of z. The 
soliton dynamics would also change to a variable speed V(z), possible due to the energy supplied by 
optical prepumping. 

This potential opportunity for controlling the pulse width, hidden in the integrable property of the 
NLS-SIT system, as explained above, was missed in earlier investigations nakazawal,nakazawa2,kakei,porsezian, 
since the initial atoms are usually taken in their ground state: \u\ 2 = 1, [P[ 2 = 0, by restricting to 
N (z) = -1. 




Figure 1: Broadening of the perturbed NLS soliton \E(z,t)\ along the fiber, moving with a constant 
speed with parameter choice k = 0.25, r] = 0.50, T = 0.28 

B. Enhanced soliton control through multiple doping 

Another promising opportunity in managing optical soliton in fiber communication, emerging also 
from the integrability of the NLS-SIT model, is overlooked completely in earlier investigations. This 
is the proposal of enhancing the effect of amplification and control of the optical soliton by replacing 
the conventional single SIT system, the only case considered in the literature, by a coupled multiple 
SIT system, using recursively the constrained integrable hierarchy in the NLS-SIT model (see Fig. 3). 
The physical meaning of coupling the NLS equation to such multi SIT system can be given through 
a novel proposal of using coupled multiple doped resonant media, in place of a single doped medium. 

For generating the governing hierarchal equations and showing their integrability, we extend Lax 
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Figure 2: a) Broadening NLS soliton pulse is controlled by a coupled NLS-SIT system with Nq(z) = 
— 0.11e rz and wq = 0.3. Variable speed of the soliton is evident from its bending in the (z,t)- plane, 
b) Additional control is achieved by coupling to a second SIT system with Nq 2 \z) = 0.11e rz , showing 
an efficient restoration of the soliton width. 

operator V(A) (I5|8h by adding more deforming terms V s um = iY^j — WoJ^Gj, linked to the M-th 
constrained hierarchy kundu09 in the NLS-SIT system, fixed at level M from the possible infinite 
sequence : j = 1,2,.... In analogy with G% (jSJ) we can express the deforming matrices Gj, through 
dipole moment pj and population inversion Nj of the j-th doped resonant medium. For explicit 
demonstration we restrict to the next higher level M = 2 in the constrained hierarchy, by considering 
only an additional SIT system to the original NLS-SIT set. Compatibility of the Lax pair thus defined 
would generate an extended set of equations given by the same deformed NLS © coupled however to 
a double SIT system 

ip lt = 2{N 1 E - w oPl - pa), iN lt = -{Ep\ - E* Pl ), 

ip2t = 2(N 2 E - w oP 2), iN 2t = ~(E P * 2 - E*p 2 ), (20) 

with induced polarization p 2 and population inversion N 2 , linked to the additional doped medium 
described by the second SIT system. We find intriguingly that, the exact soliton solution for the 
optical pulse E in this extended NLS-SIT model (|3I2Q[) . can be expressed again in the same form (fT3|) . 
where the soliton parameters are to be modified with contributions from all its interacting parts, i.e. 
from the NLS as well as from the multiple SIT system as v = v n i s +v s iti+v s it 2 , oj = uj n i s +ui s iti+u s u 2 - 
Parameters v n i s ,ui n i s and v s iti,co s itx have the same expressions as found already in (]15|16p . while the 
additional SIT contribution v s u 2 ,u} s it2, can be derived following a similar argument as f 1 1 6 [) in the form 

1 k 
v s it2 = — Im[(Xi - u>o)" 2 ]/2(» = 2—f 2 (z), 
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j 2 2 

w«*2 = i?e[(A 1 - U ;o)" 2 ]/2(^) = 2^^-/ 2 (z), (21) 

with /2(z) = - J z Nq 2 \z')<Iz' , involving an additional arbitrary function Nq 2 \z) = N2(z,t = — oo). It 
opens up therefore another novel way, hidden again in the integrable structure of the NLS-SIT system, 
for an enhanced control of the soliton width and dynamics, by adding a coupled second SIT system, 
as shown in Fig. 2b. 

This process of coupling the NLS equation to the set of multiple SIT equations can be continued 
within the framework of the integrable system, as mentioned above, creating a form of directional 
connected network with feedback, as shown in Fig. 3. In particular, as evident from the coupled 
equations (|3|20p . the input optical pulse E would influence the dipole field pj and the population 
inversion Nj in all the resonant SIT media with j = 1,2,... , M, while only p\ from the first medium 
gives feed back to the field E. On the other hand, Pj+i are coupled sequentially to pj, across the media, 
while Nj are mutually coupled only with pj from the same medium, in the multiple SIT system with 
j £ [1,M]. This network, would exhibit more and more manipulative power for control over width 
and amplification of the optical pulse, enhanced sequentially by choosing a set of initial condition 
-oo), j = 1, 2, . . . M and is based on the notion of constrained hierarchy of the 
integrable NLS-SIT system (see Fig. 3). 




Figure 3: Connected network of the NLS and the multiple SIT system with E as the input optical 
field, pj as the induced polarization and Nj as the population inversion of the j-th. doped resonant 
medium with j = 1,2, 3, . . . , M. The arrows show the directions of coupling with equations (|3l20p 
describing this network in the particular case of M = 2. Sequential enhancement of the control 
of width and amplification are predicted by this network, which is consistent with the constrained 
integrable hierarchy of the NLS-SIT system. 

This theoretical prediction, as presented schematically in Fig 3, is an experimental challenge to 
incorporate the contribution of coupled multiple SIT systems. Repeating the idea of available exper- 
imental realization of single doped fiber medium, either to a series of doped media coupled through 
induced polarization, or to multiple doping with parallel coupling in a single medium, such experi- 
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mental set up is likely to be organized. 



IV. CONCLUDING REMARKS 

Exploring the integrability of the coupled NLS-SIT system we have given novel proposals for 
controlling its solitonic pulse. The broadening problem of the optical pulse can be addressed by 
adjusting initial population inversion of the dopant atoms, linked to the soliton width, by choosing 
more general function Nq(z) > — 1 for the initial profile, in place of the traditional restriction to 
Nq(z) = — 1. This also allows amplification of the signal through initial excitation by prepumping 
energy. The controlling effect can be refined further by using another integrable property of the coupled 
NLS-SIT model given by its constrained hierarchy. The idea is to replace the conventional single SIT 
system by a network of sequentially coupled multiple SIT media with doping. Each additional SIT 
medium can bring in a new tunable function (z) > — 1, j = 2, 3, . . . in the form of initial population 
inversion of additional dopant atoms, providing more manipulative power for controlling the shape and 
dynamics of the optical soliton. One set of dopant atoms in the resonant medium is coupled to another 
set by induced polarization, with all SIT media interacting in turn with the input optical field. This 
network of interacting systems described by the constrained hierarchy of the integrable NLS-multiSIT 
equations is predicted to have enhanced control over solitonic width and amplitude, which can increase 
sequentially with the number of coupled SIT media. In such a multi-doped media requiring higher 
threshold intensity for the formation of solitonic pulse, one could possibly use a multi-level dopant like 
neodymium (Nd 3+ ), where with more than two available levels the energy can be pumped throughout 
the process, unlike in two levels, resulting to a higher gain OptCom. 

Both of our theoretical proposals with applicable potentials can be worked out analytically in 
minute details through ISM, due to the underlying integrability of the system. 
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